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Abstract 

. . , Different theoretical models which attempt to describe 
' hadronic production of heavy quarkonia are reviewed. 
. Firstly, we consider the Color Singlet Model and point 
out the large discrepancies between the theoretical pre- 
^ . dictions and the results from the Tevatron detectors. 

' Then some other models are introduced, quickly dis- 
^~~> cussed and confronted with experimental results. Fi- 
■ nally, we suggest possible ways to understand the source 
T-H , of the remaining discrepancies. 

^ ; 1 Introduction 
^ ! 

. Since the discovery of the J/tp and of some higher states 
y—i ' of charmonia, the calculation of their production rates 
has been performed. In the early 80's, there were two 
leading models for the description of the data, namely 
^ the Color Singlet ModelJCSM) H and the Color Evap- 
^ oration Model (CEM) M. Each were based on simple 
JpH, but justifiable assumptions that we shall discuss later 

Q^: oil- 

(— I Until the middle of the 90 's, the experimental results 
' were all in good agreement with these two models in 
. !^ [ any type of production but the CSM was preferred be- 
' cause of a seemingly more solid theoretical foundation. 
H ', This agreement was largely due to the fact that direct 
- - < production was not observed, so that the J/ip would re- 
sult from a variety of cascading decays. Nevertheless, 
in 1995, the CDF collaboration showed large discrepan- 
cies between the rates of ip{2S) predicted by the CSM 
and experimental results for direct production in high 
energy pp colhsions [|. The same occurred for the J/tjj 
state in 1997 when they achieved the isolation of the 
direct contribution jj]. 

This discovery gave the CEM a second life despite its 
weaker foundation compared to the CSM. Besides that, 
some other models or mechanisms were proposed to 
solve the problem. The first was the Color Octet Mech- 
anism (initially introduced for high-p^ fragmentation). 
Its key-point is that the bound state can be produced 
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in a colored state and then bleached into a singlet state 
by non-perturbative processes, these last efi^ects being 
mathematically given by non-perturbative matrix ele- 
ments 

Later, this mechanism was included in a more general 
formalism, non-relativistic QCD (NRQCD). The latter 
is based on a systematic expansion in the coupling con- 
stant and the quark velocity in the bound state, which 
for heavy quarkonia is supposed to be much less than 
the speed of light. 

The theoretical predictions based on NRQCD account 
well for all the available data from hadron colliders and 
more or less satisfactorily for data from e-p colliders {e.g. 
HERA) . The only discrepancies come from polarization 
measurements, where NRQCD predicts a transverse po- 
larization, the data clearly do not show signs of any 
polarization 

Another model was proposed by P. Hoyer and S. Peigne 
0. The basic idea is that the heavy-quark pair can 
undergo a perturbative interaction with the comoving 
color-field produced by the initial hadronic collision. In- 
troducing a new variable which parameterizes this inter- 
action, it is able to reproduce some features of the data 
which are not described by the CSM, as well as polar- 
ization measurements. 

2 The Color Singlet Model 
2.1 The model 

This model is based on several approximations or pos- 
tulates: 

• Decomposing the quarkonium production in two 
steps, first the creation of the 2 heavy quarks 
{Q & Q) and then the binding of these two quarks 
forming the meson, one postulates the factorization 
of these two processes. 

• As the scale of the first process is approximately 
M'^ +p'^ one considers it as a perturbative one. One 
can thus calculate its cross section with classical 
Feynman-diagram methods. 

• Considering only bound states of heavy quarks. 
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their velocity in the meson must be small. One 
therefore supposes that the meson is created with 
its 2 constituent quarks at rest in the meson frame. 
This is the static approximation. 

• One finally assumes that the color and the spin of 
the pair QQ don't change during the binding. As 
physical states are colorless, one requires the pair 
to be produced in a color singlet state. 

In high-energy hadronic collisions, the leading contribu- 
tion comes from a gluon fusion process. Using a meson 
production vertex function with the required tensorial 
structure and prescriptions relative to the propagators, 
one has six diagrams for the ^Si states. 





5 crossed diagrams . . . 
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In order to select the type of production, some other 
constraints can be imposed. Non-Prompt production 
(coming from h quark decay) is rejected by the detec- 
tion of a secondary vertex. The prompt but non-direct 
production (coming from radiative decay -cf. Figure 
squares and plain triangles-) is rejected by detecting the 
photon emitted during the decay. 
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Figure 2: -^B as a function of Pt for J/ijj 
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2.2 Comparison with data 

Given their quite large branching ratio into dileptons, 
the best way to detect (heavy) quarkonia is to focus on 
muon pairs and to plot their invariant mass distribution 
(cf. Figure |l] (a)&(b)). For instance, in the mass re- 
gion of the J/V', the distribution shows a maximum at 
the precise value of mj/^. Constrained fits on this dis- 
tribution for different values of kinematical parameters 
then provide us with differential cross sections relative 
to these parameters (cf. Figure ^ and || for ) . 
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Figure 1: Resonances due to J/i/i (a) and to ■ip{2S) 
(b) §. 



Figure 3: for direct T{3S) g. 

Figures |^ and |3| show the curve obtained with the CSM 
(dotted curve for Figure ||) and the measurements by 
the CDF Collaboration. In the case of the J/tp, the 
discrepancies are more or less a factor of 30, for the 
■0(25') (not shown) they reach 60 and for the T{nS) 
(only the T(3S') plot is shown in Figure |3|) the factor is 
10. 

It is therefore evident that the CSM totally fails to re- 
produce the data. The same thing happens with D0 
results. An experimental problem is thus unlikely. An- 
other important feature to note is that the electropro- 
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duction data (from HERA detectors) are up to now still 
in good agreement with CSM predictions. 



3 The Color Evaporation Model 



This model is based on the fact that agtrong > 1 for 
long range interactions or low 4- momentum transfers. 
As a consequence, the probability V that a quark pair 
undergoes many quantum fluctuations during non- 
perturbative interactions with surrounding hadronic 
matter is big. Schematically, 
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The asymptotic state (the meson) is random, still being 
colorless. Thus the probability to produce the different 
quarkonia states of a given family is supposed equal, or 
almost equal. 



Mathematically, this gives 
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where the natural value of F is the inverse of the number 
of quarkonia between the threshold 2mcM and 2mD.B- 

The new feature compared to the CSM is that the lead- 
ing contribution at low Pt consists in the following pro- 
cess 
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As already mentioned, this model is in good agreement 
with experimental data. For instance, it reproduces 
quite well the energy dependence of the cross section 
as well as its polarization. 

Nevertheless this model raises several remarks: 

• It is very phenomenological. 

• _F is in fact a free parameter, its fitted value seems 
to depend on the kinematics. 

• In order to obtain , one is tempted to introduce 
NLO contributions. This is by construction of the 
model normally included in F, and hence one may 
be double-counting. 



4 The Model of Hoyer &; Peigne. 



other related processes, e.g. open production, indeed 
suggest that there exists a comoving gluonic field pro- 
duced by the hadronic collisions. One can convince one- 
self that such field doesn't exist in general in QED, and 
in particular in leptoproduction of quarkonium (cf. Fig- 
ure ^). This could explain why the problem does not 
occur in e-p collisions as in pp collisions. Moreover, the 
re-interaction has to be perturbative because the heavy 
quark symmetry predicts that non-perturbative inter- 
actions will not change spin. In this sense the CEM 
contradicts this model and the heavy quark symmetry. 




(b) QCD (hadroproduction) 



Q 




Q 

(c) QCD (photoproduction) 

Figure 4: Due to self-interaction of gluons, the gluon 
field produced by Bremsstrahlung stays at rest in the 
rest frame of the quark pair. 

The authors of this model can thus reproduce the exper- 
imental data by fitting the variable which parameterizes 
the scattering with some assumptions on the topology 
of the field and on its polarization. 

The main success of the model is its physical content, 
its ability to explain other features than the simple pro- 
duction of quarkonium- for instance J/V' suppression 
in hadronic matter- and, as hoped, the results directly 
linked to production cross sections of vector mesons. 



5 Analysis of the theoretical un- 
certainties. 

In order to get a first idea of what could be the source 
of such discrepancies, we have undertaken an analysis 
of the theoretical uncertainties arising in the CSM. The 
first source we've considered is related to the wave func- 
tion at the origin, which enters directly the expression 
of the cross section (|l|). 

Its value is in fact extracted from the leptonic decay 
width, which writes 



This model describes the production of quarkonium 
through Hard Comover Scattering. Some features of 
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We thus find that the error on F^^ introduces an error 
of at least 10% on the cross-section. 



of relativistic wave functions and of gauge invariance in 
bound states description. 
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Then we analyzed the uncertainties due to the parton 
distribution functions (pdf). There exist two ways by 
which the pdf's can introduce errors in the cross sec- 
tion. Firstly the value of Q!s(Q^) attached with the pdf 
and secondly the pdf's themselves for the hadronic cross 
section. 

We found that the overall factor resulting from theses 
3 sources of uncertainties is about 2-3 and roughly con- 
stant for different values of Pt- 
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Figure 5: Ratio of the cross sections obtained with var- 
ious pdf's to the cross section obtained with MRS(G) 
2-94. 



6 Summary and Outlooks 

In this brief review, we saw that the CSM was unable to 
reproduce the experimental data, especially for charmo- 
nia, even if it is based on sensible approximations and 
postulates. This feature was discovered nearly ten years 
ago for the Even if NRQCD was believed to be the 
appropriate answer, according to the recent data on po- 
larization, its efficiency is now arguable. Other models, 
GEM and re-scattering model, are efhcient but we need 
more tests to reinforce their credibility. 

In this context, the evaluation of the non-static contri- 
bution of the GSM could be one of the solutions. It 
could open new paths for the understanding of this seri- 
ous problem and could help for a better understanding 
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